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Abstract 

We study a system of semilinear hyperbolic equations passively advected by 
smooth white noise in time random velocity fields. Such a system arises in mod- 
eling non-premixed isothermal turbulent flames under single-step kinetics of fuel and 
oxidizer. We derive closed equations for one-point and multi-point probability distri- 
bution functions (PDFs) and closed form analytical formulas for the one point PDF 
function, as well as the two-point PDF function under homogeneity and isotropy. 
Exact solution formulas allows us to analyze the ensemble averaged fuel/oxidizer 
concentrations and the motion of their level curves. We recover the empirical formu- 
las of combustion in the thin reaction zone limit and show that these approximate 
formulas can either underestimate or overestimate average concentrations when re- 
action zone is not tending to zero. We show that the averaged reaction rate slows 
down locally in space due to random advection induced diffusion; and that the level 
curves of ensemble averaged concentration undergo diffusion about mean locations. 
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1 Introduction 



We are interested in statistical properties of solutions of the following passively advective 
semilinear hyperbolic systems: 

u t + v(t,x,u)oV x u = f(u), (1.1) 

where u G R m , m > 1; x £ R , N > 2; / : R m — > i? m , a smooth nonlinear map; 

x, uS) = v(x) + v'(x, t, u), where v is the mean vector field (e.g. a constant vector field), 
v 1 the fluctuation random field, stationary Gaussian in space and white in time, with mean 
zero, and covariance: 

< vlix , t)^(x, t') >= D id {x - x)5(t - t'), (1.2) 

where D = (Dij)(-) is smooth. The o denotes the Stratonovitch sense of stochastic product. 
The coupling of different components of u is through the lower order term /. The specific 
forms of / in applications often permit invariant regions hence uniform maximum norm 
bounds exist on the solutions for all times. 
An example of ( |1 . 1| ) is the following system: 

Y Ft + v{t, x, uj) o V X Y F = kAY f - KY F Yo, 

Y , t + v(t, x, u) o V X Y Q = kAY - KY F Y Q , (1.3) 

where k > is diffusion constant, K > reaction rate constant, Y F (Yo) is the concentra- 
tion of fuel (oxidizer), v a prescribed random turbulent velocity field. Such a system arises 
in isothermal non-premixed turbulent combustion under single step reaction: Y+uO x — > P, 
where P is product, v stoichiometric constant, see (chapters 2 and 6), (chapters 3 
and 5), and |T7|]. We shall be concerned with the inviscid regime where k = 0. The com- 



mon quantities of interest include the ensemble averaged concentrations and reaction rate 
E[Y F ], E[Y ], E[Y F Y } (see e.g. §, section 6.7, g] chapter 5). 

It is known that these moments of solutions as well as the underlying probability distri- 
bution functions (PDF's) do not satisfy any closed equations in general. To proceed with 
an exact analytical treatment, we make the white noise in time assumption a la Kraichnan 



on the velocity field [IT], [Ml - Additionally, the sample paths of velocity field are no less 



smooth than Lipschitz continuity, see the recent work p0| , |] in this so called Batchelor 
limit. This permits us to get closed equations for PDFs for nonlinear systems of equations 
such as ( [1.1|) . To understand the analytical structures of the class of stochastic PDE prob- 
lems like ( p. . 1|) and make progress without resorting to closure approximations, the white 
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in time assumption offers a unique starting point, and helps one to gain insight. This line 
of inquiry has been pursued in many recent works jlB), 0, [0, among others. For successes 
of the white noise model in the context of the linear passive scalar model, see the in-depth 



survey by Majda and Kramer [16|. For an extension to finite time correlation from the 
white noise in time linear passive scalar model, see the recent work JI|. 

We mention in passing that similar semilinear systems also arise in bioremediation and 
transport in porous media |IS], [B3|. In [2T|, system ( |1.3| ) was derived in a fast reaction limit 
of microbial reaction when the substrate retardation factor is equal to one. However, the 



velocity v in that particular application is spatially random and independent of time |L8 



The random dependence of velocity field in time is more relevant in turbulent combustion. 
A model similar to ( |1.3|) , was recently considered by Majda and Souganidis |l7j for the 
case of two-scale periodic velocity field (oscillatory in both space and time) to understand 
the large scale effect of turbulent mixing in non-premixed combustion and condensation- 
evaporation in cloud physics. Mean field equations were obtained in the limit of large 
scale separation (homogenization limit). Here we obtain the exact equations without scale 
separation in space but with rapidly decorrelated velocity in time. 

Within the framework of the white noise model, we are able to analyze several important 
quantities ( ensemble average of solutions and reaction rate) often encountered in practice. 
We also use our exact results to understand approximate formulas in combustion literature 
on calculations of these averaged solutions. Moreover, we have made rigorous previous 
formal derivations of PDF equations by Kraichnan and coworkers flT2 |, ||], [ IU| . 



In section 2, we derive the one point PDF (P = P(u,x,t)) equation and exact solution 
formulas on averages of solutions under the initial condition P(u, x, + ) = S u=Uq m, where 
«o is a bounded nonnegative determinisitic (vector) function such that the reaction f(u) 
is integrable. In the thin reaction limit of ( |1.3D , or f(u) = except on a low dimensional 
manifold, we recover the empirical formulas in combustion (H section 6.7), and show 
how (under)overestimation can occur with these formulas when reaction zone is not thin. 
Moreover, the average reaction rate E[f(u)] decreases pointwise in time to zero with a 
much faster rate than the advection free case due to enhanced diffusion caused by random 
advection. With both a scalar model of Fisher type (f(u) = u(l — u)) and the model 
system ( |1.3| ), we show that the level curves of averaged solutions undergo diffusion about 
the mean location. 

In section 3, we derive the multipoint PDF equation and analyze its solutions. By a 
comparison argument, we show the decay of the correlation functions for large times using 
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the correlation functions of Kraichnan's passive scalar model. We also find the two-point 
PDF solution formula in closed form under isotropy and homogeneity conditions. 

Section 4 are concluding remarks. Many open issues remain for future work. One of 
them is to extend results here to the regime where velocity time correlation is a finite value 
away from zero. The other is to study more properties of level curves of concentration 
variables using PDF information. For results along this line using a different approach, see 
Constantin and Procaccia ||. 

2 One Point PDF Equation and its Solutions 
2.1 Derivation of One Point PDF Equation 

Let us first recall a representation of Gaussian random velocity v , so called proper orthog- 
onal decomposition, see []14| among others. The velocity field v(t, x, uo) is formally the time 



derivative of a cylindrical Brownian motion Wt on a separable Hilbert space H, or a Gaus- 
sian process indexed by H x R + with covariance matrix cov (W t (h), W s (h')) = min(t, s) < 
h, h' >c, h, h' G H. The underlying space H carries information of the original x variable. 
The covariance of velocity v can be regarded as a positive definite bilinear operator acting 
on H. If D = D(x, y) is in Lf oc (dx x dy), then there exists an orthogonal basis e^ n '(x) G R N 
of space H such that: 

Di,j(x,y) = £ et\x)ef\y), 1 < i, j, < N. (2.1) 

n 

Now let W^(t,co) = W t (e<ri) } which is a sequence of independent Wiener processes (with 
(jj omitted below), the representation reads: 

v(x, t,w) = Y, e (n) W dW {n \t)/dt. (2.2) 

n 

With (|2.2|) ) the original system ( |1.1|) can be put into: 



d t u = -Y J {e {n \x)-V x u)odW {n \t) + f{u)dt=-Y J g (n \x,u)odW {n \t) + f(u)dt, (2.3) 

n n 

where o refers to the Stratonovich sense of differential product. The equation ( [2.3|) with Lip- 
schitz velocity u falls into the category of stochastic flows studied in Kunita Jl3[], solutions 
exist and are unique, and can be approximated by finite dimensional stochastic ODE's. Re- 
call a well-known conversion from Stratonovich to Ito (nonanticipating) stochastic ODE's 
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(X = X(t,u) eR M , gW G R M ) 



M 



dX = f(X)dt + E 9 {m \ x ) ° dW {m \t), (2.4) 

m=l 

where W^ffis are independent Wiener processes, / and g are Lipschitz functions of X. 
Then the equivalent Ito equation is: 

i M M 

dX = {f(X) + - E g^(X) ■ V x g^\X)}dt + E g^\X)dW^\t), (2.5) 

m=l m=l 

the new drift term is the so called noise induced drift. A similar conversion for fl2.3|) 
generates the drift term: 



dyg^(y,u)V u{y) g^(x,u) = ±£ / dy {y) • V v u{y) e^(x) ■ VJ(x - y) 

^ n J ^ n J 

= Jv, • / rfy E e (n) (^) ® e (n) (y)] • vxi/W* - y) 

£ J n 

= ^ V *'/ dyD(x,y)V y u(y)6(x - y) using (2.1) 
= ^V,- (D(a:,a:)V s u(x)) 

= iv x • (D(0)V x u(a;)). (2.6) 
It follows that the Ito form of the equation is: 

cku = -E (e (n) (^) • V ( «)^ B '(i) + [/(«) + V ■ (D{0)V x u))dt, (2.7) 

n ^ 

or in PDE form: 

u t + v ■ V x u = f{u) + iv, ■ {D{0)V x u), (2.8) 

The advantage of the Ito form (|2.7|) or (|2.8|) is that the stochastic advection term has mean 
zero. 

Let ip = ip(u) : R N — > R N be a smooth nonlinear map, and Jm = Jm(u) be its Jacobian. 
Multiplying system (|1.1|) on the left by Jip, we find: 

{cp{u)) t + v o V x cp{u) = J<p{u)f{u), (2.9) 

whose equivalent Ito system of equations are with the same derivation as above: 

(p{u) t + v ■ V x ip(u) = -V x • (D(O)V^(m)) + Jcp(u)f(u). (2.10) 
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Taking ensemble mean of (2.10), we get 



E[<p(u)] t = \v x ■ (D(0)V x E[cp(u)]) + E[Jp(u)f(u)}, (2.11) 



and the last term is equal to — / (p(u)V u ■ (f(u)P(u, x, t)) du. Hence it follows from ( [2.1 1| ) 
that the P(u,x,t) satisfies the equation: 

Pt + V u • (f(u)P) = l -V x ■ (D(0)V x P). (2.12) 
We may summarize the results of this section in the following: 

Theorem 2.1 Let P(u,x,t) be the one-point probability density function of the stochastic 
solution u(x,t) of the semilinear hyperbolic system with white noise in time, spatially 
Lipschitz, incompressible, Gaussian random velocity field v(x,t) . Then P solves the closed 
equation ^2.12J . 

2.2 Applications 

2.2.1 A Model Equation with Fisher Nonlinearity 

Let us first consider a scalar model equation of the form ( |1 . 1|) with f{u) = u(l — u), the 
so called Fisher (or KPP, Kolmogorov-Petrovsky-Piskunov) nonlinearity. We shall assume 
that the advection veocity has mean equal to zero. If the advection is identically zero, then 
the solution is: 

u(x,t) = - ^ — -, 2.13 

which is between zero and one if the initial data Uq(x) is so. As t — > +oo, u — > 1 if 
x G {x : Uq(x) > 0}, and u —>■ if x G {x : Uq(x) = 0}. In particular, any characteristic 
function is invariant in time. For the random advection, the PDF equation ( |2.12| ) reads: 



P t + (n(l - u)P) u = ^D A X P, (2.14) 

with initial data P = P (u,x) = (47r<7^) -JV / 2 exp{ — [u — u (a;)) 2 /(4cr 2 )}, where u (x) is a 
bounded measurable function between zero and one, and a u is a small parameter. Such 
initial data is convenient for examining deterministic initial data by taking o u — > limit 
after we find the general formula of the PDF. 

Let us solve ( p,14| ) by taking Fourier transform in x ( justified by a suitable truncation, 
see next subsection for details) to get: 

Pt + w(l - u)P u = {2u-l- A)|£| 2 /2)A (2.15) 
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which is then integrated by the method of characteristics. The characteristic curve is: 

u ° et foi^ 
(1 - m ) + w e* 



whose inverse is: 



So: 



u + (1 — M)e t 



PoK 0(1 - % + «oe*) 2 exp{-t - ^\H} 



jo 1 — Wo + u n e l 2 



u 



,t-D \i\H/2 



P ° \u + (1 - u)e« ' V (w + (1 - w)e') 2 ' (2 ' 18) 
It follows by taking inverse Fourier transform in £ that: 

PK,^ /^""! / ^ e H^lV(-o*) (4 ^ r . /2e -(^^-«o(,)) 2 /(^)_ 
(w + (1 — w)e')" ! Vh^ 

(2.19) 

Now we calculate all moments of solution (/ = f(u) = u p , p > 1): 
t) = | duf(u)P(u,x,t) 

4yK {t , X ,y) ( du ///<"> (4w y/. e -(^«)'M, (2.20) 
J (w + (1 — w)e')^ 

where K is the heat kernel with diffusion constant Dq/2. The inner integral in ( |2.20| ) is 
equal to (via the change of variable v — u/(u+ (1 — w)e*)): 

cfc (47r ( T 2 )- Ar/2 e- { ^ Uo(y))2/(4CT ' ) /(^e7(l - v + ue*)). (2.21) 



Taking the limit a u — > 0, we find that ( |2.21|) converges to /(w (l/) e V(l ~ u o(y) + M o(?/) e *))> 
and: 

l RN -»"\->~>»/j V 1_ ' u o(2/) + w o(3/) e< 
In particular, if the initial data is a front, namely, wo(w) = x({ x : x i £ [0? oo)} = x(R+), 
fl2.22j) reduces to: 



E[f{u)] -> /_ dyK(t,x,y)f [——^L——\. (2.22) 



/ dyK(t,x,y) = (2nD t)- 1/2 / dyi e -(*i-w)7(2A>t) = ^-1/2 / e ~* 2 ^. ( 2 .23) 

jo j .a 



To probe the front motion, we look at the level set X(t) such that E[f(u)](t,X(t)) = c 6 
(0, 1). For large t, (^T23|) shows that: 



X(t) = (X 1 (t),X 2 ,---,X N ) ~ {zopD t,X 2r --,X N ), (2.24) 

where z is the unique number so that e~ z dz = The number zq > 0(< 0) if 

c > l/2(< 1/2), z Q = if c = 1/2. This implies that the random front in the average sense 
undergoes normal diffusion about its mean, in this case X\ = 0. Of course, the random level 
set {x : u(t,x) = c G (0, 1)} is more complicated and analysis of its almost sure behavior 
requires more information (multipoint statistics). 

2.2.2 The Model Combustion System 



Let us consider the one-point statistics of the two by two combustion system ( |1.3| ) with 



K — 0, K — 1, by studying the initial value problem of (|2.12| ) with f(u) = Uiii2(—l, — 1), 
and initial data: 

P ( M > X ' °) = (471-0-2)^/2 eX P{~4~2-(( ni ~ U l( X )) 2 + ( U 2 - U° 2 {X)) 2 )}, 

the smoothed delta function located at u^ix). We shall find a solution formula then 
take the limit o u — > 0. To use Fourier transform in x, we need to truncate P(u, x, 0) at 
large x, this can be done by multiplying to P(u,x,0) a smooth function ipn(x) compactly 
supported in the ball Br e R N . The truncated data is denoted by Pr(u,x, 0). Fourier 



transforming ( 2.12 ) in x gives: 

Pt ~ {u x u 2 P) ul - (u lU2 P) U2 = -A)|£| 2 P/2, (2.25) 

with initial data Pr(u, x,0), and D = D(0). Equation ( J2.25[ ) is put to the form: 

Pt ~ u x u 2 P ul - u lU2 P U2 = ( Ul +u 2 - D \Z\ 2 /2)P, (2.26) 

which we solve by the method of characteristics. The characteristic T equations: = 
—UiU 2 , i — 1,2, give solutions: 

c cie Cot c , 

U\ = -, u 2 = -, 2.27 

l + cie c °* l + Cie c o 4 v ; 



and useful relations: 



!fl = - Cl e C0 \ c = mi - u 2 , a = -^i e ( " 2 - Ul)t . (2.28) 
u 2 u 2 



7 



Along T, P obeys in view of and (|2727| ) : 



pt = («i + « 2 - a>i*i72)p = - c o - A)!^) A 



and so: 



p(e,t) = p fi («i = -^, M2 

1 + Ci 



-c A (1 + cie 



c n 2 



-( C0 +A,|S| 2 )i 



1 + d'V (l + ci) J 



(2.29) 



Using (ggp to write (pTSSQ , we find: 

(-it! + m 2 ) Mi e («2-«i)t ( M2 _ Ml ) M2 



m 2 - Mie(" 2 - ni )' ' u 2 - uxe^ U2 ~ Ux)t ' 

(«2 -U x f 



(u 2 - Mie^ 2 "" 1 )') 2 



-( Ul ~u 2 +D Q \£,\ 2 /2)t 



(2.30) 



where: 
Pr(u,0 = I dxi> R (x) 



{Anal)*/ 



2 exp{-( Ml - u\{x)) 2 /Aal - (m 2 - u° 2 (x)) 2 / 4a 2 u - % ■ x}. 

(2.31) 



Substituting fl2.31|) into (|2.30|) , and taking inverse Fourier transform, we have: 



P{u,x,t) = (2vr) 



-TV 



v dyi) R (y)- — 2-^exp{- 



(«2— Ul)t 



(-Mi + u 2 ) Mi e 
m 2 - M 1 e( U2_ ' Ul )* 



-Ay)) /Aal 



u 2 - Mie^ 2 "" 1 )' 
(u 2 - u x f Jva _ 



•Ml)* 



d£ e 



-D |£| 2 t/2-if-2/+i£-x 



(m 2 - M ie(" 2 - Ul )*) 2 

where the inner integral gives the heat kernel: (2nD t)- N / 2 exp{-|x - y\ 2 /(2D t)}. The 
exponential decay in y allows us to remove the truncation ipn by letting R — > +oo. It 
follows that: 



P(u, x, t) 



( M 2 ~ m) 2 c (n 2 - 

(m 2 - wie^ 2 "" 1 )*) 2 
i 

exp{- 



«i)t 



/ dy (2nD tr N / 2 exp{-\x - y\ 2 /(2D t)} 

JR N 



-Ml + M 2 ) Mi e 



(U2—Ul)t 



{Altai)*/* 

_ / (m 2 - Mi)M 2 
I m 2 - Mie^ 2 "" 1 )* 



m 2 - Mie^ 2- " 1 )* 

-«§(v)l /4^}. 



(2.32) 



Note that the denominator m 2 — Mi e^ 2 Ul ^ 1 being zero is not a singularity due to a similar 
term in the exponential of the Gaussian. 
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The PDF formula ( 2.32|) can be further simplified by taking o u — > for any finite time 
t. This is most conveniently done when we calculate the average quantities, E\u\ u 2 ], -^[ M i]> 
i?[u 2 ], in the limit o u — > 0. 



E[uiu 2 ] = J ui u 2 P(u, x, t) du\du 2 

dy (2nD t)- N / 2 exp{- \x - y\ 2 /(2D t)} u%) u° 2 (y) J« u° 2 , t), (2.33) 



where J = J(vi,v 2l t) = det g^'"^ , and the mapping (t>i,t> 2 ) — > (ui,u 2 ) at any time t is 
the inverse of: 

(-Mi +M 2 )Mie (u2 ~ Ul) * (U 2 -Iti)ll2 , OCM x 
Ui = / r, , v 2 = -, tt. (2.34) 

Notice that (t>i,i>2) of ( |2.34| ) satisfies the system: v itt = Viv 2 , i — 1,2, with initial data 
( v i, v 2)(fy = (wi,m 2 )- So the inverse map is the solution of the system: u i)t = —U\U 2 with 
initial data (t>i,t> 2 ), or: 

{vi - v 2 ) Vl e^- v ^ {v x - v 2 )v 2 

Ul = 7 77 , U 2 = -, r- , (2.35) 

Vie (vi-v 2 )t _ V2 Vie (vi-v 2 )t _ V2 > v / 



and the Jacobian is: 

ft ( v _ V \2 e (vi-V2)t 

J = exp{-y o (*i+*)(*)ds}= \^%_ var (2.36) 
From (|233| ) and ( gg6fl , we have: 

(2.37) 

and the total average reaction rate is: 

em* = [ ^M-^^''^ dy , (2 . 38) 

r n J r n (u\(y)e^- u >^ - u° 2 (y)) 2 1 ; 

It is natural to compare ( |2.37| ) and ( |2.38| ) with the deterministic case when say v — 0. 
Indeed, we solve two equations: u itt = —U\ u 2 , to find: 

du\{x) e dt du\{x) 
Ul u^xY* - u° 2 (xy U2 ul{x)e dt - u° 2 (x) ' 1 } 

where d = d(x) = Ui(x) — u 2 (x), u° > 0, i — 1, 2. The total reaction rate is: 

K(x) - M°(x))M(^)«2(^)e K(:r) ~"° W)t 
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which agrees exactly with ( |2.38|) ! This turns out to be true more in general on a system 
like (|1.1|) . For any smooth incompressible velocity field v, the integral J Rn g(u) dx (g a 
nonlinear function so that g(u) G L 1 (R n ) ) is independent of v by writing solution in 
characteristic variables and noting that the change of variables from x to the characteristic 
variables has unit Jacobian. However, the local (pointwise) average reaction rate E[u\U2\ 
decays faster by a factor 0(t~ N ^ 2 ) due to advection induced diffusion. 
Similarly, we find the formulas on i — 1,2: 

E[Ul] = L dv u\{y)e^-^-ul{y) 

E[U2] = L ' ^"^eK^)-^-^^ ^^^- (2 - 41) 

which are convolutions of the deterministic solutions with zero advection ( [2.39| ) with K 
the heat kernel of diffusion coefficient Do/2. If the initial data u\ = x({ x £ R '■ x i > 0}), 
u% (x) = 1 — u\(x), namely the front data, then in the large time limit, ( p. 41 ) behaves like 



and the level curves of E[ui\ = c, i = 1, 2, behave like ( |2.24| ). In other words, fronts 
undergo diffusion about the mean location X\ = in the average sense. 
To summarize, we state: 



Proposition 2.1 Consider the 2x2 model combustion system Hl.S ) obeying conditions on 
the velocity field in Theorem 1; and deterministic, nonnegative initial data (u^u^ix) G 
(L°°(R N )) 2 . Then E[ui], Efa], E[uiU2] are expressed in closed analytical form. The total 
average reaction rate E[\\uiU2\\L 1 (R N )]{t) is given by the closed form formula ( $.4Q ), the 
same as the formula when the random velocity field is absent. However, E[uiU2](x,t) is 
smaller, due to advection induced diffusion, than Uiu 2 (x,t) in the absence of random advec- 
tion. The mean concentration fronts understood in the sense of level curves of E[ui], E\a^ 
diffuse about the mean position with diffusion constant D(0). 



We remark that if the mean velocity v is present, then the PDF equation ( |2.12| ) has on the 
left hand side a term v ■ V X P. If v is a constant vector, this term can be removed by going 
to the comoving frame £ = x — vt, and the same results on reaction rates as above can be 
obtained. 



Now we draw a connection between formula ( [2. 41 ) and an empirical procedure in com- 



bustion (H, pp 247-248). Let Z = u\ — u 2 , then Z satisfy the passive scalar equation 
(without reaction), its PDF denoted by P(Z) is Gaussian with kernel K(x, y, t), in view of 
( P.12j ). If reaction occurs in a thin zone, that is, either u\ = 0, u 2 ^ or u\ ^ 0, w 2 = 0, 
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or U\ = u 2 = 0, then when Z < 0, we have Z = —u 2 , U\ = 0; when Z > 0, Z = ui, u 2 = 0. 
It follows that: 



E[ui\ = [ ZP(Z)dZ, 
Jz>o 



IZ>0 

E[u 2 ] = - [ ZP(Z) dZ. (2.42) 
Jz<o 

This provides a way to approximate the average of solutions based on the PDF function of 
the (virtual inert) passive scalar Z. 

In case of initial data: u\ = x({ x £ R N '■ x i < 0}), u 2 (x) = 1 — u1(x), we are in the 
thin domain regime. We have from ( p.41[) that: 



E[m]= [ u° 1 K(x,y,t)= [ ZP(Z)dZ, 

Ju<(>0 JZ>0 



similarly for £"[^2], hence recovering the empirical formula fl2.42j ). 

If the initial data is: u\ = x({ x ^ R N '■ x i < 0}), u 2 = e if X\ G (—5,0), u 2 = if 
x\ < —5, — 1 if X\ > 0, we have a reaction zone of thickness 5 > 0. Then E[u\] will 
differ from the thin reaction zone case above by: 



which is negative for fixed t, and small enough e. The approximation ( |2.42 ) will overesti- 
mate 

Now if u° 2 = x(i x e R N : X! > 0}), u\ = 1 if x x < 0, u\ = e, if x 1 G (0, 5), w? = 0, if 
X\ > 5, then we also have a reaction zone of thickness 5. The difference of E[u\] from the 
thin reaction zone case is: 

e(e - M°)e (e -"^ 



3*6(0,5) ee( e ~"2)* 



which is positive for small e. Hence the approximation ( ^.42 ) will underestimate E[u-\\ 



3 Multi-Point PDF Equations and Solutions 

We derive the multi-point PDF equation by imbedding n > 2 points xi, - ■ ■ ,x n G R N into 
a vector: 

X = (x ir --,x n )eR Nn , (3.1) 
and corresponding values of solutions into: 

e(X, t) = (u(x h t), • • • , u(x n , t)) G R mn , (3.2) 
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Let also: 



and: 



V(X, t) = (vfa, t), • • • , v(x n , *)), 



F(G) = (f(u( Xl , t)),---, f(u(x n , t))) E R mn . 
Then the n equations (m = u(xi,t)): 



(3.3) 
(3.4) 



Ui,t + v(xi, t, lu) ■ V x .Ui = f(Ui), 



can be written into the system: 



Q t + V ■ V x e = F(6), 



(3.5) 



which is of the same form as the original equation for u. It follows that the one point PDF 
equation of 6, or n-point PDF equation of u, is (P = P(Q,X,t)): 



Pt + Ve- ((/(ui), • • • , f(u n ))P) = 2 Vi • (D n V x P), 



(3.6) 



where D n is the n x n block matrix (D(xk — X\)) with each block an N x iV covariance 
matrix of original velocity v. 

We shall need a structure assumption on the original velocity covariance matrix D = 
D(x), x = (xi, ■ ■ ■ , x N ) e R N : 



Dij{x) = 0(0)6^ - dij{x) = D(0)«J o - - Di\x\ 



valid for the domain of small \x\, also known as the Batchelor regime, |TT 
(p.6|) can be written as: 



(3.7) 
Equation 



p t +J2v Ui -(f( Ui )p) 



i=l 



D 



~t(Y1 V ^) 2 +nE dap(xj ~ x k )d x ad : 



i=l 



3 k 



P = M n P, (3.8) 



where x® refers to the a-th component of Xj, and repeated indices mean summation. The 
operator M n is degenerate elliptic when acting on L 2 functions of x. If P is translation 
invariant, the first term on the right hand side of fl3.8|) vanishes. 

If one neglects reactions, then the last equation is of the same form as the equation 
for the Lagrangian n-particle transition probability density P^ S (X;Y), which gives the 
probability to find n Lagrangian particles at positions yp. at time s if they were at positions 
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Xk at time t < s, for k = 1, n. See ||, p. 21. There is a relation between this object and 
our multipoint probability density of the scalar amplitudes, P(Q,X, t): 

P(0; X,t) = J dYP^(X; Y)P(Q; Y, 0). (3.9) 



This relation is, in fact, just the integral solution of our equation ( |3.8|) . However, it doesn't 
hold in the general context of spatially non-Lipschitz but only Holder continuous velocities 
that was considered in ||. In fact, it holds only in the context of the "Batchelor regime" 
that we consider, where velocities are Lipschitz in space. Indeed, a simple consequence of 
our equation is that the ensemble- average "scalar energy" | / dx < 8 2 (x,t) > is conserved 
in time (without reactions), whereas this is not true in the case of non-Lipschitz velocities. 



A particular case of interest of equation (|3.8| ) holds for n = 2, if the velocity field is also 
rotationally invariant: 



Pt + V U1 • (fMP) + V U2 ■ (f(u 2 )P) = D 2 r N ~\r N+1 P r ) r . (3.10) 



Equation ( ^.10[ ) is the rigorous justification of a similar equation deduced by Kraichnan, 



(5.10) in ||12|| . He considered a single scalar with no nonlinear reaction. His equation was 
written for the probability density Q(A,k,t) of the 2-point difference A = u — v of the 
scalar at wavenumber k. Kraichnan's proposed equation read 

Q t = D 2 [{k{kQ k ) k + NkQ k ] (3.11) 

in the special case where molecular diffusion vanishes. We obtain this equation by changing 
variables in formula ( |3.12j ) (written for the case of m = 1 and / == 0) from u and v to 
A and U = (u + v)/2, by integrating out the second variable U, and by then Fourier 
transforming from space variable r to wavenumber k. 

For the system ( |1.3j ), if we consider rotationally invariant (isotropic) solutions for the 
two point PDF P, then ( p. 10 ) becomes: 



P t + {-u x u 2 P) Ul + (- Ul u 2 P) U2 + (-v lV2 P) Vl + {-v x v 2 P) V2 = D 2 r N -\r N+1 P r ) r , (3.12) 

where D 2 is a positive constant; we have used (161,1*2) to denote the solutions at the first 
point, (vi,v 2 ) at the second point, and their separation distance is r. Assuming that the 
initial data and solutions decay rapidly at space infinities, we multiply (|3.8|) when n = 2 
by uiu 2 viv 2 , integrate over (u,v), to get (E 2 = E[uiu 2 viv 2 ] for short): 

(E 2 ) t = D 2 r N ~ 1 (r N+1 E 2 ) r - J u\u 2 v\v 2 {u\ + u 2 + v 1 + v 2 ) P(u,v, x,t) dudv, (3.13) 
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implying that < E 2 < G = G(r, t), where G is the solution of: 

G t = D 2 r 1 - N (r N+1 G r ) r , (3.14) 

which is the equation for two point correlation function of Kraichnan's passive scalar model 
in the free decay regime, see [P0| |, and 0. 

To solve either (|3. 14 ) or ( 3.12| ), let us make the change of variables: £ = iVlogr, then 
the operator: 

D 2 r 1 - N {r N+1 d r -) r = D 2 iV 2 (% + 9 e ) • . (3.15) 
Solution to (|3TD is then (D 3 = D 2 N 2 ): 

G = (AnD 3 t)-^ 2 J dr ] ex P {-^- V + D 3 t) 2 /(AD 3 t)}G (e^ N ), 

= {AnDzty 1 ' 2 ex V {-{N\og{r/s) + D 3 tf/{AD 3 t)}G Q {s)Ns~ l ds, (3.16) 
Jo 

which shows exponential decay rate on any compact set of r away from zero for fast decaying 
Go(s). G(0,t) is however conserved due to smoothness of G near r = 0, or not decaying 
to zero within a neighborhood of 0(cq®^ e~ D ' il / N ) to leading order, cq a constant > 1. 
The behavior of G implies that f( XjV -\ e n x n E[ui(x)u 2 (x)ux(y)u 2 (y)] dxdy converges to zero 
exponentially as t — > oo, Q e R a compact domain. Using L 2 decay property of the 



semigroup generated by the operator M n [14"], the decay of all correlator functions follow 
without isotropy. 



We solve (|3.12|) in (u, v, £, t) by combining Fourier transform in £ and method of char- 
acteristics in (u, v, t) as before. The result is: 

(u 2 - M X ) 2 (v 2 - v t ) 2 



P(u,v,U) 



- J (u 2 -u 1 +v 2 -vi)t 



(u 2 - Ml e(" 2 - U1 )*) 2 (v 2 - vie( V2 - V1 ^) 2 

[ dy (AitD 3 t)- 1 / 2 exp{-|£ -y + D 3 t\ 2 /(W 3 t)} 

JR 1 



■P, 



'(-ui + m 2 ) Mi e { " 2 - Ul) * ( 



U 2 — U\)U 2 



)u 2 {-v l + v 2 )v 1 e^-^ t ( 



''2 - e y/ N 



\ u 2 - u 1 e { - U2 ~ Ul "i t ' u 2 - Mxe^ 2 -" 1 )* ' v 2 - v 1 e { - V2 - v ^ t ' v 2 - v 1 e( V2 ~ Vl ') t '' ) 
where Pq = Po(u, v, r) is the initial PDF data. Back to the original variable r, we have the 
two-point PDF formula: 



[u 2 - U\, 



(v 2 - Vif 



^(u2—Ul+V2 — V-i)t 



V ' r ' ^ \u 2 - Ul e^- U ^) 2 (v 2 - v^-v^ty- 

/ Ns- 1 ds(47iD 3 t)- 1/2 exp{-\N\og(r/s)+D 3 t\ 2 /(4D 3 t)} 
Jo 



Po 



-ui + u 2 ) u\ e*" 2 Ul )* (u 2 — Ui)u 2 (—vi+v 2 )vie lyV ' 2 (v 2 — Vi)u 2 



u 2 — U\e 



(U2-Ui)t 



u 2 — U\e 



{Ul-Ux)t ' y 2 _ Vie ('02-Vi)t ' y 2 _ VlC (V2-Vl)t ' 



(3.17) 
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Summarizing the results of this section, we state: 

Theorem 3.1 Let P(u±, u n , x, t) be the n-point probability density function of the solu- 
tion u(x,t) of system with advection velocity satisfying the conditions of Theorem 1. 



Then P obeys the closed equation ( \3.o\ ). If the separations rij = \xi — Xj\ are all small for 



i,j = 1, ...,n, then the equation takes the special form (\3.8{). If in addition v is rotationally 



invariant, the 2-point PDF P(ui,U2, x,t) obeys the simplified equation \3.1ty and the closed 
form analytical solution for P is given by ( \3. 1 % ) for the case of deterministic initial data 
of the system (fJT; 



4 Concluding Remarks 

Under the white noise assumption of incompressible advection velocity field, PDF equations 
are closed for the solution of the semilinear hyperbolic system ( |1 . 1|) . Exact formulas for one- 
point and multi-point PDFs (under isotropy) on solutions of the isothermal non-premixed 
turbulent flame system ( |1.3| ) can be derived. From exact formulas, one is able to recover 
empirical formulas in combustion on averaged solutions in the thin reaction zone limit, 
and analyze the effect of finite domain zones on the empirical formulas. It is found that 
they can either overestimate or underestimate averaged solutions. It is also found that 
fronts understood as level curves of averaged concentration variables, diffuse about the 
mean position in the ensemble average sense. The front diffusion about mean location 
is known in the distribution sense in the context of the Burgers' equation and other one 



dimensional convex scalar conservation laws, see fl22| . It is conceivable that this picture 
is valid for turbulent combustion fronts when the white noise assumption is relaxed. A 
future study with assistance of numerical computation appears feasible. It will also be 
interesting to study more geometric properties of random level curves using higher order 
PDFs, such as its roughness and dimensions in the large time limit for front data in several 
space dimensions. See for results in this direction with a different approach. 
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